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We analyze electronically excited nuclear wave functions and their coherence when subjecting 
a molecule to the action of natural, pulsed incoherent solar-like light, and to that of ultrashort 
coherent light assumed to have the same center frequencies and spectral bandwidths. Specifically, we 
compute the spatio-temporal dependence of the excited wave packets and their electronic coherence 
for these two types of light sources, on different electronic potential energy surfaces. The resultant 
excited state wave functions are shown to be qualitatively different, reflecting the light source from 
which they originated. In addition, electronic coherence is found to decay significantly faster for 
incoherent light than for coherent ultrafast excitation, for both continuum and bound wave packets. 
These results confirm that the dynamics observed in studies using ultrashort coherent pulses are 
not relevant to naturally occurring solar-induced processes such as photosynthesis and vision. 



Introduction 

An increasing number of studies in recent years have 
identified surprisingly long-lived coherences in biologi- 
cal systems, especially in photosynthesic light harvesting 
molecules [THS] . Such long-lived coherences were detected 
in studies such as two-colour photon echo spectroscopy 
[3J 3], angle-resolved coherent optical wave-mixing [3], 
and phase-stabilized 2D electronic spectroscopy [7]. Fol- 
lowing these findings, it has been suggested that quan- 
tum coherence may play an important role in biological 
processes, for example in the proton and energy trans- 
fer that follow the photo-excitation step. A challenging 
question is to what extent the dynamics observed in such 
experiments, which utilize coherent laser light sources, 
are relevant for processes induced by incoherent sources 
such as sunlight, even when the two types of light share 
the same center frequencies and spectral bandwidths. 

Several previous studies have analyzed various effects 
of incoherent excitation sources on the resultant molec- 
ular dynamics. These include an analysis of the excited 
state survival probability as a function of the incoher- 
ence of the light [8], effects of incoherent fields on the 
photoisomerization yields [9] , a quantum-optical formula- 
tion of the state of the molecule after photon-absorption 
[TU] . studies of open system dynamics relevant to pho- 
tosynthetic complexes [TT], and general features of the 
response of open systems to incoherent excitation [12]. 
However, none of these studies looked directly at the spa- 
tial shape of the wave packets and the associated tempo- 
ral evolution and coherence properties resulting from ex- 
citation with incoherent vs. coherent light sources. Such 
studies are particularly important in comparing the na- 
ture of the excited quantum state prepared by these light 
sources. 

In this paper we concentrate on the actual spatio- 
temporal shape of the wave functions resulting from 
photo-excitation, a primary entity that effectively brings 



out the differences between the two types of light-induced 
excitations mentioned above. In addition, we consider 
the photo-excitation of a model molecule into a super- 
position of two electronic states, and contrast the cal- 
culated coherences induced by incoherent solar-like light 
with the outcome of excitation by an ultrashort coher- 
ent pulse with the same center frequency and spectral 
bandwidth. 

Our main theoretical tool is a time-dependent 
quantum-theoretic calculation of the different nuclear 
wave packets resulting from the photo-excitation step. 
Optical incoherence is simulated by introducing random 
jumps in the phases and center frequencies of the exci- 
tation pulse, with the degree of optical incoherence be- 
ing controlled by the frequency and amplitudes of these 
random jumps. Both the spatial structure and time de- 
pendence of the resultant wave packet dynamics is then 
analyzed for a variety of potential energy surfaces (PES). 



Theory of "random" and coherent wave packet 
dynamics 

Consider first the case in which a nuclear wave packet 
is photo-excited to a single excited PES. The molecule 
which is assumed to be initially in the ground electronic 
state is excited by a light source, treated as a classi- 
cal time- varying electric field, to a dissociative or bound 
excited state. The time-evolution of the total molecu- 
lar quantum state can be described by a (stochastic) 
Schrodinger's equation, written in the adiabatic (Born 
Oppenheimer) electronic basis as 

ifi^¥(R,t)=g(t)*(R,t), (1) 

where 
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with "fg (^ x ) being the nuclear wave packet in the ground 
(excited) electronic manifold. The total Hamiltonian ma- 
trix consists of a molecular part and an electric-dipole 
Matter-Radiation interaction part 



H(i) = H M + H MR (t), 
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The assumption that the molecule is initially in the 
ground electronic manifold, ]J/(R, to) = M^s^'^-M 5 re- 
duces Eq. Q to 

<MRt) = k{t m T) U x (t-T)D(R)E(T)U g (T-t )^ g (R,to), 

(8) 

where 



(2) 



The nuclear kinetic energy operator Tn and the ground 
and excited PES functions W g>x of the above depend on 
the set of nuclear coordinates R. The Matter-Radiation 
interaction term is given as, 

V MR (R,t) = - J r x *(R iqe )d(R,q e )^ g (R,q e )dq e -eE(t), 

(3) 

where q e is the set of electronic coordinates, E(t) = 
£h{t) cos(wt+$) is the electric field, d(R, q e ) is the dipole 
moment, and e is the unit polarization vector. 

It is most convenient to solve the Schrodinger equation 
in the interaction representation (with the R variable 
temporarily suppressed for brevity) 



1 for f > t 
h(t - t) = { \ for t = r 
for f < t 



(9) 



is the Heaviside function. The physical picture that 
emerges is that no population is excited before the ar- 
rival of the delta pulse, and that exactly at t — r, where 
h(0) = 1/2, and U e (r — t) = 1, we have that 



ih 



dt 
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U f (t - t ) • *(to) 
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where to represents some initial time long before the onset 
of E MR (t), and 

gf R (t) = g(t-t )E MR (t)]i(t-t a ) , 

with the (interaction representation) evolution operators 
defined as 





U x (t f -U) 



m - ti ) = e-m M (tf-t*) = (U a (tf- t i) 

(5) 

Assuming the laser field intensity is weak enough so 
that first-order perturbation theory is valid, the state in 
the interaction representation assumes the form 

U t (t-to)-*(t) = *(to) + 4 tuf^-moW, (6) 



or in the Schrodinger's representation, 



*(t) = U(t - < ) 



- /*Hf fl (t)dt' 



*(*o) 



(7) 



In order to gain physical insight we now consider ex- 
citation by a E{T)5{t - r) pulse, for which V MR (R, t) = 
D(R)E(T)5(t - r), where 



£>(R) 
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Thus, at t — t we form on the excited PES a "replica" of 
the ground state wave packet times D(R)E(t) /2ih, that 
evolves at t > r as 



* x (R,f >t) = U x (t-T)V X (R,T). 



(11) 



Using the above insight we now expand any time- 
dependent function for the electric field as a series of 
delta pulses, 



E(t) = / E{r)6{t - r)dr. 



(12) 



and generalize Eq. (18]) for a general pulse to obtain that 



* B (R,i) 



* C^(t-r) 



ih 



D(R)E(T)U g (T-t )V g (R, t )dr. 

(13) 

The excited wave packet is seen to be a sum of replicas 
launched by all the delta pulses up till time t, evolving 
according to the excited state Hamiltonian H x [T5J [TS] . 
Written in differential form 



ih-* x {R,t) = H x * x (R,t)+D{R)E(i)U g (t-t )y g (R,t ), 

(14) 

which is an inhomogeneous Schrodinger's equation for 
the excited state nuclear wave function ^ x (R, t) whose 
source term is the field times the dipole operator times 
the ground electronic state. 



The great advantage of Eq. ( 14 ) is that it allows us 



to easily realize varying degrees of incoherences in the 
light sources by introducing random jumps in the cen- 
ter frequency and the carrier-envelope phase of a select 



3 



CD 



_CD 
LU 



3 

CD 



O 
_CD 

LU 




100 150 
Time (fs) 



200 250 




20 30 
Time (fs) 




12 



15 



Time (fs) 



0.05 0.1 0.15 

Energy (Hartree) 



FIG. 1: Upper left panel - a sample realization of incoherent light. Upper right panel - magnified view. The center wave- 
length is ~455 nm (corresponding to center energy of 0.1 Hartree). The pulse has envelope function £h{t) ~ [sin(-7rf/200 fs)] 01 
for < t < 200 fs, and is zero elsewhere. The incoherent light simulates solar radiation by introducing, every 7 fs on the 
average, a random phase jump in the [— tt, it] range, and a center frequency shift in the ±0.0175 Hartrees range. Lower left 
panel - A coherent pulse whose envelope is £t(i) ~ exp \—(t — 6 fs) 2 /2(1.6 fs) 2 ] . Lower right panel - The frequency profile 
of the pulses of coherent and incoherent light, demonstrating that both light sources share the same center frequency and the 
same spectral bandwidth. 



number of the 5 function components of the electric field. 
Specifically, we write the magnitude of each 8 function 
component of Eq. (fl2l) as 



(15) 



E(t)=£ l (t) cos{[w L + Aw(t)]t + $(t)}. 



The real-valued functions Aw(r) and $(t), which are 
taken to be constant (zero) in time for coherent fields, 
serve as random variables for incoherent light. 



Results and Discussion 

Consider then the results of numerical simulations by 
pulsed coherent and incoherent light. Figure [TJ displays 
the coherent pulse, and one specific realization of the 
incoherent light pulse, in which we have chosen similar 
average values for the center frequency and bandwidth, 
with rate of jumps of the center frequency shifts and 
phase interruptions modelling typical solar values. The 
total energy flux, J_ \E(t)\ 2 dt, is kept the same for the 
two types of pulses. 

As a model molecular system we consider the PES of 
the molecular ion. We examined the photo-excitation 



from the ground la g electronic state, to two (dissocia- 
tive or bound) excited states of the la u symmetry. We 
assumed that these two states are decoupled from one 
another in order to reduce the dynamics to that of a sum 
of the individual ground to excited state transitions of 
Eq. (14). Although this is a model system, the general 



results that emerge are expected to apply to a wide va- 
riety of systems. 

We assume that at to = the ground state wave func- 
tion is a ro- vibrational energy eigenstate with zero angu- 
lar momentum. The transition dipole function for a lin- 
early polarized excitation field along the z-axis, is given 
as, -D(R) = D(R)Yi(6) where R is the internuclear sep- 
aration. We model \D(R) \ R/2 a.u. as in [T3] , 

Writing the J = 1 wave function as a product of vi- 
brational and angular parts ^(R, t) = Y®(6, <j})i/j x (R, t), 
and using atomic units, transforms Eq. (|l"3|) to 



ip x (R,t) = I Ui =1 {t-T,R)D{R)E{T)e- lE * T Tp g (R,Q)dT. 

(16) 

The excited state propagator for the J — 1 channel is 
given as 



XJt=\t -t,R) = exp [-iH J x =\R){t - t) 



(17) 
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FIG. 2: The spatio-temporal plots of the modulus-squared radial wave functions of the vibrational wave packets, u x ' (R,t). 
Top row is for purely repulsive excited state-1 and state-2 potentials, and bottom row is for bound Morse wells. The functional 
forms for both types of PES correspond to the first and third situation in Figure [3] 
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(18) 



with fi being the reduced mass associated with the R 
coordinate. 



Using Eq. (14) and Eq. (16), we have 
dip x (R,t) 



dt 



-iHi- 1 {R)tP x {R,t)+D{R)E{t)e~ lE o t i) g {R). 

(19) 

Finally transforming to radial wave functions 
Ug : x(R,t) = Ripg tX (R,t), we obtain the radial ana- 



log of the inhomogeneous Schrodinger equation, 



d 

dt 



u x (R,t) = -i 



-Id 2 
2[i dR 2 
D(R)E(t)e~ 



+ V X J =\R) 



u x {R,t) 



iE„t 



Ug{R). 



(20) 



Modelling the electric field E(t) as in Fig. [I] and choos- 
ing u g (R) as the v = 5 vibrational state, we used Eq. 



( 20 1 to propagate the pair of excited radial wave func- 
tions u x 1 \R, t) and u x 2 \R,t) from t = to 300 fs, us- 
ing time-steps of 0.003 fs and radial-steps of 0.02 Bohrs. 
We analyze the motion of the u x ^{R,t) and u x (R, t) 
pair that has originated from the same ground electronic 
state for two cases: 1) Continuum dynamics for a pair 
of purely repulsive exponentially-decaying excited poten- 
tials. 2) Bound state dynamics for a pair of Morse po- 
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tentials, where a heavier reduced mass, p, = 5 x 918 a.u. 
is used. Sample plots of the wave packets, for these two 
cases, subject to the coherent and incoherent modes of 
excitation, are presented in Figs. [2j In the case of the 
incoherent excitation, the results of a single random re- 
alization is shown. 

The completed wave packets resulting from incoher- 
ent light excitation are seen to be choppy and highly 
unstructured. In the case of the repulsive potentials, fol- 
lowing the excitation step, these wave packets spread out 
in configuration space much more rapidly, while covering 
much larger extensions, than do the coherently excited 
cases. The situation is even more dramatically seen in 
the bound state excitation cases where the wave pack- 
ets resulting from the incoherent light excitation exhibit 
a significantly smaller number of coherent oscillations as 
compared to excitation by coherent pulses, and a highly 
irregularly structured wave function is seen. 

The above results show characteristics of the excited 
state wave functions on individual electronic states. Also 
of interest, in particular with respect to experimental 
studies of long-lived coherences, is the persistence of elec- 
tronic coherence, a property of several electronic levels. 
To examine this we calculate the electronic coherence 
given by pi,2(i), the off-diagonal element of the density 
matrix of the two electronic states, traced over the nu- 
clear spatial coordinate, 



Pl,2(t) 



(R,t)u£>(R,t)dR. 



(21) 



Because the coherence is obtained via first-order pertur- 
bation theory, it is convenient to normalize pi,2{t) by 
dividing it by pi.i(t) + p2,2(t), the density matrix trace, 



giving 



C(t) 



pi At) 



+ 02,2 (*) 



(22) 



In this way we factor out the effects of the total excitation 
yield that are of no interest for this comparative study. 
Note also that Eq. (21 1, by averaging over vibrations, 



includes the effect of decoherence, due to the vibrational 
degree of freedom, on the electronic coherences. The re- 
sultant normalized coherences C(t) are plotted in Fig. 
[3j for different choices of the PES for the two electronic 
states, under the action of the two types of fields. 

These plots clearly show that for incoherent field ex- 
citations the magnitude of the coherence C (t) decays at 
a much faster rate than that of the coherent pulse, ap- 
proaching a much lower asymptotic value than that of 
the coherent case. This is especially so, as shown in the 
middle row of Fig. [3j when coherence C(t) is not signif- 
icant, to begin with, in the coherent case, due to small 
Franck-Condon overlap regions of the two surfaces. In 
the case of the bound Morse potential wells, clear and 
gradual decoherence of the two wave packets is observed 



for the case of coherent excitation. This is in strong con- 
trast with the incoherent pulse excitation, where the co- 
herence C{t) decays almost immediately. 

In either of the three cases the final electronic coher- 
ence induced by the incoherent field is a small fraction 
of the coherence generated by the coherent case. Thus 
we find that the incoherence of the light eliminates "long 
lived coherences" when the excitation is carried out with 
pulses of solar-like incoherent sources. 

Several comments are in order. First, note that deco- 
herence observed here does not include additional envi- 
ronmental effects, since the molecule is isolated, i.e. the 
system is closed. Rather, the observed fast coherence loss 
arises , in the incoherent light case, from the nature of the 
two wave packets as they move away from one another. 
That is, the molecular vibration serves as the decohering 
environment [16j . and it is particularly effective in this 
case due to the highly unstructured and choppy nuclear 
wave functions created by incoherent light. If, in addi- 
tion, coupling to an external environment is present, the 
decay of electronic coherence will also reflect this cou- 
pling |12j . As a consequence, long lived coherences are 
not expected to persist in realistic open systems irradi- 
ated with incoherent sources. 

Second, note that the incoherent light source used 
above acts over a 200 fs time scale. Thus, although 
the 7 fs coherence time of solar radiation is represented 
in this computation, there are two significant differences 
between the results of this computation, and that which 
would result from solar irradiation, which acts over far 
longer times, and that is is effectively CW. First, any 
pulse possesses some degree of coherence due to the pulse 
envelope [5]. Hence, the "incoherent pulse" used here 
possesses more coherence than would natural solar radi- 
ation that is incident for minutes or longer. Second, at 
such long times, as discussed elsewhere 0110], long time 
excitation of isolated molecules using natural incoherent 
light leads, when averaged over realizations, to station- 
ary eigenstates of the Hamiltonian that do not evolve in 
time. Relaxation in open systems also leads to mixtures 
of stationary states |12j . 



Conclusions 

This paper has shown that the nature of the light- 
induced wave function, as well as the rates of electronic 
coherence in molecules, strongly reflect the coherence 
properties of the incident radiation. Coherent light pulses 
produce well localized wave packets whereas incoherent 
pulses produce irregularly structured coordinate space 
densities. In addition, we have calculated electronic co- 
herences in the photo-excitation of a molecule into a su- 
perposition of two electronic states, using the two types 
of electromagnetic pulses with different degrees of coher- 
ence. The results show decreased electronic coherence 
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FIG. 3: Left column - PES and the v = 5 vibrational wave function. Top panel Repulsive PES of the form, Wx\R) = 
exp[-{R- R^)/a (i) ] +w£>, i = 1,2; where W (1) = 0.1 Hartree, R ( 1} = 3.3 Bohr, a (1) = 1.0 Bohr, w£ } = 0; W (2) = 0.1 
Hartree, R^ = 3.0 Bohr, a (2) = 1.2 Bohr, w£ ] = 0.005 Hartree. Middle panel: W,* 1 ' = 0.1 Hartree, R^ = 3.3 Bohr, 
a (1) = 1.0 Bohr, w£> = Hartree; W (2) = 0.1 Hartree, Rf> = 2.5 Bohr, a (2) = 1.0 Bohr, wj? = 0.01 Hartree. Bottom 
panel: W# '(R) = Wtf } [l - e -(fl-4* > )/ [l<i) ] 2 + w£\ i = 1,2; where W (1) = 0.1 Hartree, R^ = 5.0 Bohr, a (1) = 2.5 Bohr, 
W£ ] = Hartree; W (2) = 0.1 Hartree, R^ = 5.01 Bohr, a (2) = 2.53 Bohr, wH ] = 0.01 Hartree. Right column - The 
absolute value of p(t), the normalized electronic coherence resulting from coherent and incoherent light excitation for the three 
different PES of the left column. The results for the incoherent light were averaged over a set of 10 random-jump realizations. 



values and much faster decoherence rates when the ex- 
citation is conducted with pulsed incoherent light than 
with excitation with coherent pulses of the same center 
frequency and spectral bandwidth. Thus, even in the 
absence of external environmental effects, the random 
character of the incoherent light enhances decoherence. 
This adds support to the conclusion [10l[T2l[T7] that long- 
lived electronic coherences observed in biological systems 
are a result of the use of coherent ultrashort pulses for 
molecular excitation in these experiments and that such 
coherences should not be expected in natural processes 



induced by solar radiation. One would anticipate that 
the differences noted here between coherent and incoher- 
ent excitation would be observable in two dimensional 
photon echo experiments performed with coherent vs. in- 
coherent light |18j . 
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